function[MM]=G_matrix_mod2D(xout,yout,Nout,xb,yb,N,mu,epsilon)
for i=1:Nout
    xi=[xout(i), yout(i)];
    [M2N]=G_matrix_2XN(xi,xb,yb,N,mu,epsilon);
    MM((i-1)*2+1,:)=M2N(1,:);
    MM((i-1)*2+2,:)=M2N(2,:);
end